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We characterize the class of all physical operations that transform Gaussian states to Gaussian 
states. We show that this class coincides with that of all operations which can be performed on 
Gaussian states using linear optical elements and homodyne measurements. For bipartite systems 
we characterize the processes which can be implemented by local operations and classical communi- 
cation, as well as those that can be implemented using positive partial transpose preserving maps. 
As an application, we show that Gaussian states cannot be distilled by local Gaussian operations 
and classical communication. We also define and characterize positive (but not completely positive) 
Gaussian maps. 
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I. INTRODUCTION 



Many applications in the field of quantum informa- 
tion require the ability of preparing general states and of 
performing arbitrary transformations with them. How- 
ever, there are physical systems where the set of states 
that can be generated as well as the transformations that 
can be implemented are very restricted. For example, 
in quantum optical systems linear transformations in- 
volving beam splitters, phase plates, homodyne measure- 
ments and polarizers are readily implemented whereas 
more general ones can only be performed with a low ef- 
ficiency. Moreover, with these tools and a squeezer one 
can generate only a small class of states, the so-called 
Gaussian states. Despite this fact, in these systems a 
surprising richness of quantum information protocols has 
been found within the realm of linear optics: entangle- 
ment generation [Q , teleportation j| , key distribution || , 
quantum error correction cloning some of which 
have already been implemented [^| . For the moment it is 
not known how (or whether) important operations such 
as entanglement distillation or a (useful) depolarization 
of continuous variable states can be implemented with 
linear optics § [§ |, [h| || || . This raises the question 
which transformations can in general be realized by the 
concatenation of such operations, i.e. the tools currently 
available in the lab. 

A partial answer to this question is contained in a 
mathematical paper written in the seventies Jl~3|] . There 
a subclass of operations that transform Gaussian states 
into Gaussian states have been mathematically charac- 
terized, namely those which are trace-preserving. The 
relation of these maps to some of the experimentally fea- 
sible operations has been discussed by Eisert and Plenio 
in |l2| , and used by these authors to derive a criterion for 
the interconvertibility of two-mode Gaussian states under 
feasible local transformations. Unfortunately, the set of 
operations considered in Ref. jl3) does not include mea- 
surements (which are not trace-preserving). However, 
one of the main strengths of linear optics is the highly ef- 
ficient measurement of the quadrature observables X and 
P which homodyne detection affords. Furthermore, mea- 



surements followed by classical communication have been 
seen to be an essential ingredient in certain basic quan- 
tum information protocols such as quantum teleportation 
Jlij or entanglement distillation [T^). Thus it is impor- 
tant to find a mathematical formulation of the physical 
actions that can be applied to Gaussian states in which 
also not trace-preserving operations (measurements) are 
included |6|. 

In this work we give a full answer to this problem by 
providing a simple description of all operations of this 
sort. As an application we discuss the question of dis- 
tillation of Gaussian states by local Gaussian operations 
and classical communication (LOGCC). In this context, 
it was recently shown (TtJ that the entanglement of a 
symmetric two-mode Gaussian state of two parties can- 
not be increased with the help of another copy of that 
state and a homodyne measurement. Here we give a 
proof that distillation is not possible for an arbitrary 
number of modes per site, general Gaussian states and 
general Gaussian operations. 

This paper is organized as follows. In Section |lj we 
fix our notation for Gaussian states, setting the stage for 
the results given in the following sections. Section |n] 
contains the main results of this paper: we give a charac- 
terization of all completely positive maps that transform 
Gaussian states into Gaussian states (Gaussian opera- 
tions). We show that they can all be implemented with 
the currently available means. We derive a simple, com- 
pact form of these maps. In Section IV we consider Gaus- 
sian operations on bipartite systems; we classify them 
with respect to their locality and separability properties. 
As an application of the methods introduced before, we 
show that Gaussian states cannot be distilled by using 
Gaussian LOCC. In Section |y| we describe positive Gaus- 
sian maps and characterize them completely. Appendix 
^ contains some material on a new entanglement mea- 



sure for Gaussian states introduced and used in Sec. IV 



II. GAUSSIAN STATES 

We consider the Hilbert space of n harmonic oscilla- 
tors H = Li^lR"). A Gaussian state is described by 



2 



a density operator p whose characteristic function |lq] 
X P (x) ■= tv[pW(x)] is a Gaussian function in x € IPv 
(or, equivalently, by a Gaussian Wigner function, which 
is related to \ by a symplectic Fourier transformation). 
The operators 

W(x) = exp [~ix T R] , (1) 

are the Weyl operators (displacement operators) and 
R = {X 1 ,P 1 ,X 2 ,...,P n ), with [X k ,Pi] = i8 kl . Since 
the density operator is bounded, we can write without 
loss of generality Jl9| 

p = TT"" / dxe~i xT ' yx+tdTx W(x), (2) 

where we have used tr[W(x)] = ir n 5(x) Jll| to normalize 
p. Occasionally, we will denote p as in Eq. (^) by p 7>d . 
The matrix 7 = j T > i J n is a 2n x 2n real matrix called 
correlation matrix (CM) and d is a 2n real vector called 
displacement. These two quantities fully characterize the 
Gaussian state p. The symplectic matrix J n is 

fe=l ^ ' 

(We will omit the index whenever there is no risk of con- 
fusion.) Note that the CM usually contains all the inter- 
esting information about the properties of the state which 
are useful for quantum information in general; in partic- 
ular, the entanglement properties of a Gaussian state is 
solely determined by its CM. Thus, in some parts of this 
paper in order to simplify the notation we will omit the 
displacement from our discussions when it does not play 
a relevant role. 

Of course, Gaussian states can be also defined for com- 
posite systems, e.g. those whose density operators are in 
B(H)®B(H). An important example of a Gaussian state 
is the maximally entangled state $ p0[ . The state $ is 
the limit r — » 00 of Gaussian states (n identical two- 
mode squeezed states) with CM 

- ( % % ) < 4 > 

where A r — coshrl and C r — sinhrA are 2n x 2n square 
matrices, and 

A = diag(l,-l,l,-l,...,-l). (5) 

The density operator of $ is a projector on the improper 
state vector |3>) 12 oc J2k>o 1^)2- 

III. GAUSSIAN OPERATIONS 

Physical actions are mathematically characterized in 
terms of completely positive (cp) maps acting on the cor- 
responding density operators. The best way of charac- 
terizing them is using the isomorphism between cp maps 
(physical actions) and positive operators (unnormalized 
states) 



A. General form of Gaussian operations 

Given a cp map £ acting on bounded operators B{TL), 
we define the positive operator E G B(H) <8> B(H) as 
follows 

£i2 = (£®l)(|S>ia(*|)- (6) 

This equation has a direct physical meaning. It tells us 
that given £, we can always obtain the state E by prepar- 
ing a maximally entangled state and acting with the map 
on the second subsystem. Conversely, given the state 
E E B(H) <8> B(H) and a state p £ B(H), if we measure 
the second subsystem of E and p in the "Bell basis" , i.e., 
an orthonormal basis of maximally entangled states con- 
taining |$), and obtain the result corresponding to the 
state 1$), then the resulting state is £ (p). Thus, given 
the state E we can always implement (probabilistically) 
the map £ provided we can perform Bell measurements. 
This can be viewed as "teleporting p through the gate £" 
pgj . In formulas, we have 

£(p) cx tr 2 [Eflp 2 ] = tr 23 (E 12P3 |$) 23 ($|) (7) 

Thus, given a cp map we can generate the corresponding 
state and given the state we can physically implement £ 
(probabilistically) (23). 

Now we define a Gaussian completely positive (g-cp) 
map, Q by the properties that both Q and 1 <g> Q map 
Gaussian states to Gaussian states. With the help of the 
isomorphism it is straight forward to characterize them. 
First, we use the fact that the state |4>) appearing in (||) 
is Gaussian, so that the corresponding operator G must 
be Gaussian |24| ]. We write it as 

G= [ dxe-^ xTTx+tDTx - c W{x). (8) 

Clearly G > iff T, D, C are real and Y > iJ; that is, Q 
is a g-cp map if the Gaussian operator isomorphic to Q is 
described by a proper CM T and, conversely, to each such 
operator corresponds a g-cp map Q. Since all Gaussian 
states can be generated (e.g., from the vacuum state) 
by unitary Gaussian operations and discarding subsys- 
tems |2f| this shows that all Gaussian operations can 
be implemented by these means plus Bell measurements 
(homodyne detection). 

Now, we determine the action of the map Q on a gen- 
eral Gaussian state p in terms of G. Apart from normal- 
ization, we have 

G ■■ P~f,d^ P-y'.d', (9) 
and we find for 7', d'\ 

i = f 1 -f 12?r ?— f£, (10a) 
1 2 + 7 

d! = D 1 +f 12 ^—{D 2 + d), (10b) 
Y 2 + 7 
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where we have denoted 



1 12 



Tl2 

r 2 



D 



Di 
D 2 



and 



f = (i ® A)r(i © A). 



(11) 



(12) 



Thus, we have that all g-cp maps on 13(H) are character- 
ized by a correlation matrix Y > iJi n and a displacement 
vector D. 

In order to derive Eqs. we just use Eq. ((?]), re- 

placing E\2 by G and /?2 by a Gaussian state p 7i <; as in 
Eq. (|2|). In evaluating the trace we use the commuta- 
tion relation W{x)W{y) = e^ 2xTjy W(x + y) p) and 
tr[W(x)] <x S(x) to obtain 



(13) 



py t d> oc tr 2 (G T2 p 7!( i), 
with 7',d' as in Eqs. ©. 

B. Examples 



How to interpret the operation described by (T,D)7 
To better understand what actions (r, D) describe, we 
now briefly discuss how the familiar Gaussian operations 
are contained in our formalism. To do this, we apply 
these operations to the first subsystem of the maximally 
entangled Gaussian state <& with CM T — lim r ^oo j(r), 
with 7(r) as in Eq. (||). 

Obviously, the identity operation corresponds to the 
maximally entangled state i.e., to Y = lim r 7(r),_D = 
0. Now, performing a displacement operation on $ 
leaves the CM unchanged, but produces a displacement 
D = (Di,Q). Now we turn to the trace-preserving 
g-cp maps considered in [O. These describe all ac- 
tions that can be performed on p by first adding an- 
cillary systems in Gaussian states, then performing uni- 
tary Gaussian transformations on the whole system, and 
finally discarding the ancillas. On the level of CMs 
these operations were shown to be described by 7 1— > 
M T jM + N. A Weyl operator W{x) is mapped to 
(det M)- l e- l l ixT ( M ~ 1 *> TNM ~ lx W(M~ 1 x) by these oper- 
ations. It then follows that the Gaussian operator that 
corresponds to this operation has the CM 



r = lim 



M T A r M + N M T C r 
C r M A r 



Using formulas (0) gives 7' = lim r M T A r M + N - 
M T C r A(AA r A + 7)- 1 AC r Af. For r -> 00 we have 
(AA r A + 7)" 1 -> A; 1 ~K X lK X + o(coshr)- 3 which 
yields the desired result M. 

Finally, we consider an example of Gaussian measure- 
ments. The typical measurement is homodyne detec- 
tion, which realizes the von Neumann measurement of 
the operator X. It has been shown before Wf\ that 



with the use of an ancillary system and a beam split- 
ter, homodyne measurements may be used to realize the 
generalized measurement corresponding to the positive- 
operator-valued measure (POVM) {|a)(a|,a € <D}, 
where a) is a coherent state, i.e. in the language of CMs 
a state with CM 7=1 and displacement d = (Rea, Ima). 
Since every other pure Gaussian state can be obtained 
from I a) by Gaussian unitaries, this implies that all 
POVMs of the form {\j,d){j,d\ : d G R 2 ™} can be 
performed with homodyne detection and suitable pre- 
processing. 

To see which CM T corresponds to the measurement of 
7, d) (7, d\ we apply it to <E>. In order to get an interesting 
result, we perform only a partial measurement, i.e. we 
consider an (n + m)-mode system and measure only the 
last m modes. The corresponding operator is 

£7=<7,d|$|7,d>, 

where I7, d) describes a pure Gaussian state of m modes. 
Expressing $ in terms of Weyl operators, replacing 
(7, d\ W(x) I7, d) by e~ x / ix i x ~ ld x and integrating over 
the modes measured we obtain a CM 



(14) 



and a displacement D = (0,0,rf) T . Note that the first 
row corresponds to system "1", while the second and 
third row refer to system "2" . This represents a straight 
forward generalization of the situation considered in Sec- 
tion [II to maps which decrease the number of modes 

(0) for a n + to mode Gaussian 




present. Evaluating Eq 
state with CM 7 



A C 
C T B 



we obtain fit 



i = A-C 



d> = lc' 

2 B + lp 



(15a) 



(15b) 



which corresponds to the change in CM derived in jlTj for 
projections into pure Gaussian states. Homodyne detec- 
tion itself represents the limiting case in which the CM 
7 P becomes infinitely squeezed. In this limit, the inverse 
in Eqs. ( |l5|) is to be understood as the pseudo-inverse 
(inverse on the range). 

In general, noise-free Gaussian operations (unitaries 
and von Neumann measurements) correspond to pure 
state CMs T and noise added to the CM describing the 
operation directly translates into noise added to the out- 
put state, i.e., we have that r' = Y + P > T implies 
that Gr'(7) > Gr(7)V7 > iJ. To see this, consider the 
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operation Qy-> and write G in Eq. ( |l3| ) as a mixture of 
states with CM T. This shows that the state Gr'{p~t,d) is 
a Gaussian mixture of states Gr(Py,d+d 2 ) displaced by di, 
where x — (di,^) are distributed according to a prob- 
ability distribution proportional to cxp(— l/Ax T P~ 1 x). 
But since a displacement cfe of the input state does only 
affect the displacement of the output state [cf. Eqs. ([To|)], 
it follows that Qr> (pj,d) is nothing but a Gaussian mix- 
ture of states Gr(Pi,d) displaced by some value y which 
is distributed according to a Gaussian distribution with 
covariance depending on P and 7. Thus the operation 
Qt' could be realized (for known 7) by first performing 
Qy and then performing random displacements to add 
the appropriate noise, which proves the assertion. Since 
displacements can be done locally, this becomes particu- 
larly useful in the discussion of entanglement distillation 
with Gaussian means below. 



C. Deterministic Operations 

In general, the transformation Eq. (j^) is not trace- 
preserving. This is related to the fact that we have con- 
sidered only one of the possible Bell measurements in 
Eq. (0). The projector |$) («&| can be extended to a 
POVM by considering all displacements W(x) |<I>), x £ 
R 4 ™. Using the second relation in Eq. ([?]) it is easy 
to see that if GtAp) = tr 23 (-Bi2P3 1$} 2 3 ($1) triCn 
GtAp) = ^2 3 [E 12 p 3 W{-D) |$) 23 ($| W(D)], i.e. D can 
be understood as the (continuous) output of the Bell- 
measurement implementing Q. 

Note that D has no influence on the CM of the result- 
ing state. Hence, provided 7 and L are known, Q can be 
turned into a trace-preserving operation by postprocess- 
ing: conditional on the measurement result {D\,D-i) the 
corresponding displacement -Di+ri2(r2+7) _1 -C>2 can be 
undone, leading to a deterministic transformation that 
maps every Gaussian state p 7i d to p 1 , f 12 (f 2 +7)- 1 ci w ith 
certainty. Note that this is true even if 7 is a state on 
a multipartite system, since displacements can be done 
locally. It is a curious feature of Gaussian operations 
that even measurements do not change the CM non- 
deterministically. 



IV. BIPARTITE SYSTEMS. APPLICATIONS 

In this section we consider Gaussian maps Q on bipar- 
tite systems. In this situation it is interesting to distin- 
guish whether Q can be implemented with local opera- 
tions on the subsystems A and B (possibly enhanced by 
classical communication, LOCC) or whether interaction 
between the system is necessary. Our formalism yields a 
very convenient form for any local Gaussian operations, 
and allows to determine the nonlocal properties for any 
given Gaussian map. 

As an application we use our formalism to show that 
entanglement distillation is not possible with Gaussian 



means. This extends the results of [jl7| to any number 
of modes, all kinds of Gaussian operations, and all kinds 
of Gaussian states. Note that, there are other means 
of performing distillation which do not require Gaussian 
maps as long as a Kerr non-linearity or photodetection 
are available |t], [|. However, at the moment these proto- 
cols still pose considerable experimental challenges and 
none has been implemented to date. 



Local Gaussian maps assisted by classical 
communication 



To determine whether Q can be implemented with lo- 
cal operations on the subsystems A and B (possibly en- 
hanced by classical communication, LOCC) or whether 
interaction between the system is necessary we can use 
the the ideas of Ref. (2j| (which extend the Jamiolkowski 
isomorphism to bipartite systems). This allows to read 
off the answer to this question from the Gaussian state G 
isomorphic to Q. If G is separable, then it can be gener- 
ated by local action and classical communication. Note 
that following the discussion at the end of the previous 
section this implementation can be done deterministically 
provided the CM of the state on which we act is known. 

If G is entangled, two cases can be distinguished: is G 
has positive partial transpose (ppt), the corresponding 
map can be implemented with a so-called ppt-preserving 
channel, otherwise full-fledged quantum interaction be- 
tween A and B is needed. 

The separability criterion for Gaussian states |3(J 
allows us to decide for every given map, whether it is sep- 
arable 1 30 , ppt-preserving P9| or neither. Moreover, the 



characterization of separable CMs given in 1 3 1 , namely 
that 7 is separable iff 37^4, jb > iJ such that 7 > ja®1b 
implies that - except for added correlated noise - all 
Gaussian LOCC operations are of product form. 



B. Gaussian states cannot be distilled with 
Gaussian local operations and classical 
communication 

Entanglement distillation is a process in which two sep- 
arate parties A and B transform a large number of copies 
of a bipartite mixed entangled state pab (jointly written 

as P®b) a state ip^AW which, as n goes to infinity 
approaches a pure maximally entangled state. To this 
end, A and B are allowed to perform arbitrary local op- 
erations (correlated by classical communication) on their 
respective part of p'ab- In the following we show that 
such a process is not possible when pab is Gaussian and 
only Gaussian operations are allowed. 

We consider a bipartite system composed of subsys- 
tems A and B and a partially entangled Gaussian state 
with CM jab and want to check if a separable Gaus- 
sian map can increase the entanglement. To this end 
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we define a simple function V(j) to quantify the entan- 
glement of a general bipartite CM. Let V(7ab) be the 
largest value p < 1 such that jab > p{lA © 7b), for 
some CMs ja and 7b . Note that for a maximally entan- 
gled Gaussian state V = [this and further properties 
of V(j) which are used in the following are proved in 
Appendix [X| , and therefore the goal of a Gaussian distil- 
lation protocol would be to decrease V(jab)- In general 
one would allow an arbitrary number copies of the state 
with CM 7ab i i- e - a state with the CM 0^ =1 7ab- As 
shown in App. |A|, V does not change when adding more 
copies of the same state. The question then is, whether 
there is a local Gaussian operation Q that produces from 
these states an output state j' A b wu -h V{i'ab) < V(iab) 
or even allows to reach V — > in the limit of infinitely 
many copies. 

In the following we show that with local Gaussian op- 
erations it is impossible to decrease V at all. Our proof 
makes no assumptions on the size or type of the entan- 
gled state considered or the local Gaussian operations 
performed. In particular, it covers any number of copies 
of a n x n-mode input state. 

We consider a separable Gaussian completely positive 
map acting on two systems A and B. As discussed above, 
the action of such a map on the correlation matrix ^ab 
is completely characterized by another correlation ma- 
trix r acting on an extended space of systems A, A', B 
and B'. The fact that the map is separable implies that 
r = r aa' © ^ b b' + P, where P is a positive matrix . In 
light of the discussion at the end of Subsec. 



[II B 



this 

means that T can be implemented by first performing 
the (completely uncorrelated) operation corresponding to 
Taa> (BTbb> and then performing (classically correlated) 
random displacements of the resulting state according to 
a probability distribution depending on P and the CM 
jab of the input state. Since these displacements do not 
increase the entanglement, we can concentrate on the ef- 
fect of the product transformation Taa' © Tbb'- 
Now let 



AA' 





C a \ 




(B 1 


C b \ 


Wa 


A 2 J 


Tbb' — 




B 2 ) 



(16) 



Let us denote by Jab ^ ne correlation matrix of A and B 
after the action of the map T aa 1 © ^b b' ■ Then we have 
Iab > Ra® Rb, where 



R A = A 1 -C A -~- 



Ao 



-Cl>A 1 -C A ^ 



1 



PI A 



A 2 +piJ 



AT 

(17) 

and similarly for Rb- Now, we use that Taa' > 0, i(J®J) 
since T is a CM and therefore Yaa' > p[iJ © (— iJ)]- 
This implies that the RHS of Eq. ([l7|) is > ipJ and we 
immediately obtain that j' A = 1/pRa is a correlation 
matrix. From this follows that j' A b — p(i'a © 7b) an< ^ 
therefore V {j' AB ) > V (jab )■ 

What does this imply for distillation? First, it proves 
that the maximally entangled state cannot be approached 



even asymptotically (i.e., in the limit when initially in- 
finitely many copies of 7ab are available). This fol- 
lows directly from the fact (cf. App. [A]) that ^(7ab © 
Iab) — V(jab), i-e., V is invariant when adding more 
copies of the same resource. So entanglement distilla- 
tion of Gaussian states with Gaussian means is impossi- 
ble. More generally, V(jab) puts a bound on all state 
transformations that can be achieved by Gaussian LOCC 
[and even Gaussian LOCC supplemented by an unlim- 
ited amount of auxiliary entangled Gaussian states of 
V( Paux ) > V( 7AB )}. 

However, the result still leaves room for interesting en- 
tanglement transformations with Gaussian means. The 
best thing that could happen - respecting the bound 
set by Vmin — V(jab) - is to have the pure entangled 
state with V(p pure ) = V m in- Our proof does not rule 
out the possibility of "entanglement purification", i.e., of 
transforming a large number of Gaussian mixed entan- 
gled states into a (asymptotically) pure entangled state 
with the same value of V with Gaussian means. First 
calculations indicate that this might indeed be possible. 
These results will be reported elsewhere. 



V. GAUSSIAN POSITIVE MAPS 

In this Section we show how to extend the approach 
presented in Section [II to include Gaussian positive but 

To this end, 



not completely positive (g-p) maps. 10 tms end, we 
first define the set of Gaussian operators, generalizing 
Gaussian density matrices to not self-adjoint operators. 
Every operator A G BCH) is completely determined by 
Xa{x) ■— ti:[AW(x)] |18(]. It follows that A may be writ- 
ten in terms of xa as |19[ 



A = 7T-" / cLxxa(x)W(-x) 1 
Wc define the set of Gaussian operators on TL by 



(18) 



Q(H) := {A : A = [ dxe 



yx+ib x—c 



W(x)}, (19) 



where j T = 7 6 M 2 „(C),6 G (D 2 ™, and c G C. It is 
straight forward to check that A is (1) bounded iff Re7 > 
0; (2) self-adjoint iff I1117 = 0,Im& = 0,lmc = 0; and 
(3) positive iff s.a. and 7 > i J; and (4) tr(A) — Tr n e~ c , 
where tr[H / (x)] = Tr n S(x),x G K 2 ™ was used. To prove 
(3), consider one mode, 7 = gt 2 - The corresponding 
operator is \/g > diagonal in the number basis (for 
g > 1 it describes the well-known thermal states of a 
field mode) and the eigenvalues (n\ A g \n) are seen to be 
all positive iff g > 1, otherwise odd numbers correspond 
to negative eigenvalues. Finally, recall that all selfadjoint 
Gaussian operators can be transformed into (a tensor 
product of Gaussian operators of) that form by quasifrce 
unitaries, performing the normal mode decomposition of 
p or, equivalently, the symplectic diagonalization of 7, 
which concludes the proof. 
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Now we turn to linear maps on B(H.). Generalizing 
Section III we define a Gaussian map as a linear map 
Q : B(H) — > B(H') that maps Gaussian operators to 
Gaussian operators, i.e., Q\Q(H)] C Q{TL'). Again, we 
can use the isomorphism of |GM to show that all Gaus- 
sian maps on B{TL) correspond to Gaussian operators on 
B(H) ® B(7i). i.e., they may be described by a matrix T, 
a vector D, and a phase/normalization constant C. Then 
we can use Eq. (0) to calculate how the Gaussian map 
Q corresponding to (T, D, C) acts on the Weyl operator 
W(x). One finds 

g[W{x)] = f dye-i& T r(l)+' lDT (l)~ c W(y), 



where T, D belong to the partial transpose of the operator 
G isomorphic to Q. 

One quickly convinces oneself that the map Q is self- 
adjoint, that is Q(A>y = Q{A) iff G is selfadjoint, i.e., iff 
Imr = 0, ImD = 0, ImC = 0. 

A Gaussian map Q is called a Gaussian positive map 
(g-p map) iff it maps Q(H)+ to Q(H)+. In terms of 
matrices this means that Q is g-positive iff Q (7) > iJVj > 
iJ. Expressing the action of G through its matrix T this 
is equivalent to the condition 

f 1 - Tu— — f 12 > iJ V7 > U. 

r 2 + 7 

Inverted this inequality it is seen to be equivalent to I\ > 
i J and 

7 + f 2 - f £ f 12 > V7 > iJ, (20) 

1 1 — iJ 



M 



which can now be written in a 7-independent way as 

{z\M + iJ)z, z\M - iJ)z] > 0. (21) 



To see that this is equivalent to Cond. fl20| ) note that for 
any z = z r + izi G C 2 ™ there exists a symplectic map S 
such that z^SS T z = \iz^Jz\ — 2\zJ'jz i \ = c 2 ; this can 
be seen immediately by extending {si = z r /c, s 2 = Zi/c} 
to a symplectic basis {sk} and defining S by Ssk = £k, 
where refers to the canonical basis. Therefore, if cond. 
( pp| ) is to hold for all 7, it holds in particular for 7 = S T S. 

The minimum of these maxima can for given T be ef- 
ficiently sought numerically, thus providing a practical 
characterization of all positive Gaussian maps. We em- 
phasize, that such a practical characterization of positive 
maps is not currently available for general maps on d- 
level systems. 

Condition @ says that the matrix T of a positive 
map can be such that neither N :— T — iJ > nor N := 
r — iJ > 0, but there may be no vector x G (D 2ra such that 
both x T Nx and x T Nx are negative. Examples are (i) cp 



maps (r > iJ, cf. Section III); (ii) maps for which T 
iJ but r > iJ. These are decomposable positive maps 



(such as transposition); (iii) Gaussian maps for which 
neither N nor N is positive can also be constructed; these 
and the case of non- decomposable g-p maps, i.e. those 
that are not derived from transposition and the relation 
of g-positivity to the usual notion of positivity will be 
discussed elsewhere [||). 



VI. CONCLUSIONS 

We have characterized all the physical actions that can 
be performed using linear optics, squeezers, and homo- 
dyne measurements. We have also characterized those 
that can be implemented with LOCC and those that can 
be implemented using ppt-preserving maps. We have 
used the methods developed in the preceding sections 
to show that Gaussian states cannot be distilled by local 
Gaussian operations and classical communication. 

Finally we have extended the definitions given before 
to general linear maps that map Gaussian states to Gaus- 
sian states and provided a complete characterization of 
positive Gaussian maps. This emphasizes that Gaussian 
states are worth studying not only because of their ex- 
perimental relevance (which will reduce as non-Gaussian 
states become more accessible), but also on mathemati- 
cal grounds that this class of states is simple enough to 
derive strong results while being large enough to encom- 
pass most (if not all) aspects of entanglement. 



APPENDIX A: PROPERTIES OF V(j) 

In this section we collect a number of useful properties 
of the quantity V{^) introduced in Sec. [Tv|p. V is defined 
for bipartite CMs 7 (or, equivalently, for Gaussian states 
P~f,d) as 



Vh) := 



V 7A,7 



{p < 1 : J > p{jA ® Jb)} ■ (Al) 



(1) V for more than one state: 

V(7©7') = rnin{V(7),V-(7')}. 



(A2) 



that is, V does not decrease when several entangled states 
are joined together. Rather, V of the combined state is 
given by the smallest V of the individual states. 

To see this, let v = V(j),v' = V(j'). Clearly, 
V(j 7') > min{w, v'} since by definition of V we have 
7 © 7' > v(j A © 7b) © v'(i A © i B ) > vam{v, v'}{"i A © 
1b ®l' A ©7b)- On the other hand y(7©7') <min{i>,i/} 
since 7 © 7' > V(j © 7')(7aa' © Jbb 1 ) also holds for the 
reduced states with subsystems AB or A'B' traced out. 
More generally, it follows that V(®klk) — min{V(7fc)}. 



(2) An upper bound for V: 

V(p 1 ) < min{A mi „(7), 1}, 



(A3) 



where Xmm{l) is the smallest symplectic eigenvalue 
(smaller than 1) of the CM 7 of the partially transposed 



7 



state p^ A [^3). \ m in < 1 is necessary and sufficient for 
the corresponding state to have a non-positive partial 
transpose. 

(3) V for the maximally entangled state |<£>): 

V($) = lim VfaM) = 0, (A4) 

r— »oo 

since A mm (7(r)) = e~ r . 



(4) V and negativity [34|: For 1 x N systems (i.e., in 
systems where no ppt-entanglement exists we also 
have V(p) > \ m in since in that case j/X m in has positive 
partial transpose and therefore is separable. This 
shows that V(p) is related to the negativity measure of 
entanglement [f54|, and for 1 x N systems — log 2 [V(p 7 )] 
coincides with the log-negativity (up to a factor). In 
contrast, for ppt-entangled Gaussian states (3^] V{^) 
is strictly smaller than 1, while the negativity of such 
states is zero (X m in > !)■ 

(5) V as a Gaussian measure of entanglement: 

We have seen that V(j) does not decrease under local 
Gaussian operations. Hence it can be considered a 



measure of entanglement for Gaussian states. In view 
of (4), we see that V does quantify both npt and ppt 
Gaussian entanglement - in contrast to most other 
measures of entanglement calculated to date. 

(6) V is computable: It is worth pointing out that 
V("f) is also computable, as one can use the separability 
criterion derived in p0| to find the largest p for which 
7/p is separable. 

Note added: Upon completion of this work we learned 
that Jaromfr Fiurasek J35| | independently arrived at a 
similar description of general Gaussian operations and, 
in particular, of Gaussian LOCCs. 
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